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Abstract 

The exact solution of a particular form of the stationary state generalized Fokker-Planck equa- 
tions, which is given under certain conditions by the classical Tsallis distribution, is compared with 
the solution of the MAXENT equations obtained using the classical Tsallis entropy. The solutions 
only agree provided the Tsallis parameter, q, is no longer taken to be constant. 
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The nonextensive entropic measure proposed by Tsallis [j], Q] introduces a parameter, 
q, which is not defined 0, 0]. It has been argued that perhaps this parameter should be 
constant, if not universally, at least for classes of dynamical systems 5H7| . Attempts have 
been made to set limits on the value of this parameter For the practitioners this has 
been accepted de facto and calculations involving the Tsallis entropy have generally used one 
piece of data to determine the value of q which is then fixed over the range of interest j^, fl . 
However, for the generalized Fokker-Planck(GFP) equation j^], dissipative optical lattices jlo| 
and more recently in the analysis of relativistic heavy ion collision 11] q is most definitely 
not a constant, unless as in the latter case, complicated scaling arguments are invoked. 

In the present letter we consider the GFP equations 12J where an exact solution (both 
static and time dependent) has been found and is shown under certain circumstances to be 
equivalent to the Tsallis classical distribution functions. In the static case we formulate and 
solve the equivalent MAXENT problem[l3| for different values of the equation of constraint. 
We solve for q in a manner suggested earlier [l^] and find the numerical dependence of q on 
the Lagrange multiplier associated with the equation of constraint. 

Consider the following non-linear one dimensional GFP equation 



9F d rw 1 d 2 [F 2 -i] 



(1) 



where F is the distribution function, Q is the diffusion coefficient, and K(x) is the drift 
coefficient which determines the potential: 



V{x) 



K(x)dx. 
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The particular power q — 2 is chosen in accordance with the quite general discussion of the 
generalized Bogulubov inequalities 15|] which points out that systems which obey Tsallis 
statistics exhibit abrupt changes at q — 2. 

Assuming the boundary conditions F — > for x — > ±oo, the differential equation govern- 
ing stationary solutions = 0) can be written as 



1 d(F 2 ~ q " 



(3) 



By performing the redefinition u = F 2 q , separating variables and integrating, one obtains 
the following stationary solution 



F(x) = D[l-P(l-q)V{x)]r- 



(4) 



where D is a positive (undetermined) integration constant. Eq. fl4]) is precisely the form of 
the distribution function obtained from Tsallis Statistics [l, 



9] . One can verify by substitution 



that Eq.fll]) is a solution to Eq. (3) provided 
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Clearly, /3 is a function of q. It should be noted that this distribution need not be normalized 
to 1. Also, because the GFP Equation is nonlinear, for a particular value of /3 and q there 
is only one value of D that solves Eq. fl3]). Therefore, D is also an (undetermined) function 
of /3 and q. Note that for F(x) > 0, /3 > and q< 2. In the limit q — > 1 the GFP reduces 
to the standard Fokker-Planck equation, f3 — > ^ and the stationary state adopts the usual 
form 

F(x) = Dexp[-(2/Q)V(x) (6) 
The stationary state distribution can also be obtained from the classical Tsallis entropy 

S q = -^(F-F q ) (7) 

via the MAXENT equation Q 

5 F S q = (8) 

along with the equation of constraint 

Tr[F q V] = C (9) 

The solution of these equations yields the classical Tsallis distribution given in Eq.flU) for 
any allowed value of q. In order to simultaneously determine q an additional equation 

dS, 
dq 



1/3 



(10) 



must be solvedjl^]. Note the distribution obtained is not normalized. 

For the simple case where V is a constant, the phase space integrals in the trace each 
produce a constant factor which can be absorbed into C, yielding 

F q V = C (11) 

By solving Eq. ( TTTj) for (3 and substituting it into Eq. (8), we numerically solve Eq. ( TTUj) 
for q. FIG. ([B is a plot of vs. q for V = 1. 
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FIG. 1: vs. q for V=l, and = ^D^ 1 

This curve can be fitted by the function (3 = where D Q = 1.96 and Q=l. As 

already mentioned, D in Eq. (JSJ) is an undetermined function of /3, and therefore q. To 

i 

match the MaxEnt and GFP solutions, one must identify D = D Q {2 — g) 9-1 . 

Since q is a variable one cannot choose its value from one data point and assume it 
remains fixed over the entire range of interest. To further illustrate this point consider C(/3) 
as given by Eq. ( ITT]) , for constant V. For a thermodynamic system this can be thought of 
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FIG. 2: C vs. p for V=l, /3(g) = ^Df 1 



as the thermal response of the system. FIG. (j2J) plots C vs. (3 for fixed values of q as well 
as for (3 = qD^ -1 . Clearly the results obtained both in magnitude and shape do not agree 
well except in the region around the value of /3 corresponding fixed chosen value of q. 

Our results indicate that some care should be taken when the Maxent equations are 
used with the Tsallis entropic measure to describe experimental data. Although we have 
only considered the stationary solutions of a particular form of the one dimensional GFP 
equation, the MAXENT solution is only consistent with the exact solution when q is no 
longer treated as a constant. Similar behavior is also to be expected in the solutions of the 
time dependent GFP equation. Even if a scaling function can be found for every dynamical 
situation, the functional dependence of the unsealed Tsallis parameter should be determined 
as a function of any Lagrange mutipliers associated with the equations of constraint. 
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